Phase Transitions Between Topologically Distinct Gapped Phases in Isotropic Spin 

Ladders 
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We consider various two-leg ladder models exhibiting gapped phases. All of these phases have 
short-ranged valence bond ground states, and they all exhibit string order. However, we show 
that short-ranged valence bond ground states divide into two topologically distinct classes, and 
as a consequence, there exist two topologically distinct types of string order. Therefore, not all 
gapped phases belong to the same universality class. We show that phase transitions occur when 
we interpolate between models belonging to different topological classes, and we study the nature 
of these transitions. 



I. INTRODUCTION 

Even though the spectrum|-of the spin-1/2 Heisenberg 
chain was obtained by Bethall almost seventy years ago 
using his famous Ansatz, low-dimensional spin systems 
are still a strong area of activity, full of surprises and 
puzzles. Agmajor source of this activity was Haldane's 
conjecturep which predicted that isotropic antiferromag- 
netic Heisenberg chains with integer spin have a gapped 
spectrum, while chains with half-integer spin have a gap- 
less spectrum. There has been considerable theoretical 
and experimental evidence in support of Haldane's con- 
jecture; it is probably appropriate to call it a theorem, 
despite the lack of a rigorous mathematical proof.u In- 
cidentally, the gapped phase in integer spin Heisenberg 
chains has come to be known as the Haldane phase. 

A breakthrough in understanding the nature of the 
Haldane phase came when it was realized that one can 
go without a phase transition from the spin-1 Heisenberg 
chain to the AKLT modeljj where the ground state is 
made up solely of nearest neighbor valence bonds. The 
Haldane gap is thus related to the energy needed to break 
short-ranged valence bonds. 
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FIG. 1. (a) Typical configuration of the spin-1 chain, (b) 
Neel order after removing all sites with Sf = 0. 

Another important step was when den Nijs and Rom- 
melse identified a, hidden order in the Haldane phase of 
the spin-1 chain.B They showed that although the sites 
with Sf = 1, 0, —1 are not well ordered in position, their 
sequence is ordered in the way shown schematically in 
Fig. ^. That is, if we remove all sites with Sf — 0, the 
remaining sites have Neel order. The order parameter 
which reveals this hidden order is the non-local string 



order parameter 



O" == - lim ( Sf exp 



U J2 sn sf\ , (1) 



where S°' is the spin-1 operator at site i, and a — x,y, z. 

A further impetus for the study of low-dimensional 
spin systems wasrffiven recently by the discovery of spin- 
ladder materials.cl Since the spin-1/2 Heisenberg chain 
has a gapless excitation spectrum with spin-spin correla- 
tion functions exhibiting power-law behavior, it initially 
came as a surprise when it was found that the two-leg lad- 
der had a gapped spectrum with exponentially decaying 
spin-spin correlations,!] while the gapless spectrum sur- 
vived in the three-leg ladder. Thus ladders could have a 
gapped or gapless spectrum, depending on the number of 
legs. More specifically, Heisenberg ladders with an even 
number of legs have a gapped spectrum, while ladders 
with an odd number of legs have a gapless spectrum. 

The appearance of a gapped spectrum for even-legged 
ladders and a gapless spectrum for odd-legged ladders 
is highly reminiscent of Haldane's conjecture for spin 
chains. Therefore, a natural question arises: Is the 
gapped phase in spin ladders related to the Haldane 
phase in spin chains? In particular, does the gapped 
phase exhibit string order? For the rest of the paper, we 
will restrict our discussion to the two-leg ladder. 

For the case of ferromagnetic interchain coupling, it 
is clear that the two-leg ladder can be equivalent to the 
spin-1 chain, with the two spins on a rung forming an 
effective S = 1. Such a phase is directly related to the 
Haldane phase, and the system has string order. Simi- 
larly, if the interchain coupling is antiferromagnetic and 
along plaquette diagonals, when the interchain coupling 
is equal to the coupling along the chains, the model is 
in fact the composite spin representation of the spin-1 
chain; thus, the gapped phase is equivalent to the Hal- 
dane phase. (See Sec. II for a more detailed discussion 
of the composite spin representation.) Besides the gap, 
these models exhibit another characteristic feature of the 
Haldane phase. Namely, the ground state is unique if pe- 
riodic boundary conditions are used, but it is four-fold 
degenerate for open boundary conditions. 



A gap appears in the excitation spectrum of spin lad- 
ders for other types of interchain couphng as well, namely 
for antiferromagnetic coupling along the rungs or ferro- 
magnetic coupling along plaquette diagonals. Whether 
these gapped phases are related to the Haldane phase 
is much less obvious. For example, for antiferromagnetic 
interchain coupling along the rungs, the zeroth order pic- 
ture is given by studying the limit in which the inter- 
chain coupling is much larger than the coupling along 
the chains (i.e., J± ^ J). In this limit, the ground state 
is well described by a product of rung singlets, with an 
energy gap ~ J± to break a singlet and form a triplet ex- 
citation. Here, the ground state is unique, irrespective of 
whether open or periodic boundary conditions are used. 
Nevertheless, as was demostrated by WhitcjS the antifer- 
romagnetic ladder can be transformed continuously to a 
model (seemingly) equivalent to the composite spin rep- 
resentation of the spin-1 chain by switching on an irrel- 
evant further neighbor coupling. Consequently, the anti- 
ferromagnetic ladder is also related, an. some way, to the 
Haldane phase and has string order .ErB 

There are also other types of ladders exhibiting spin 
gapped phases, so the question of string order in ladders 
and the relationship of the gapped phases to the Haldane 
phase becomes even more interesting. In particular, the 
spin-1/2 chain with second neighbor coupling is often rep- 
resented as a two-leg zig-zag ladder.Ej For a particular 
value |-©f the couplings, the so-called Majumdar-Ghosh 
point JIj the ground state is known exactly. It is dou- 
bly degenerate and the excitation spectrum is gapped; 
each ground state consists of a sequence of independent 
singlets. However, it was shown that the-^ajumdar- 
Ghosh groupd state has perfect string order .113 It has also 
been showr£3 that the Majumdar-Ghosh model can be 
smoothly connected to the ladder with strong ferromag- 
netic rung coupling, which is equivalent to a spin-1 chain, 
without a phase transition. 

Since the ground state of the Majumdar-Ghosh model 
consists of decoupled singlets, similar to the ground state 
of the ladder with strong antiferromagnetic coupling 
along the rungs, one might get the impression that all 
spin-gapped phases can be smoothly connected to each 
other. By this we mean that one can go from one model 
to another by continuously varying the model parameters 
without undergoing a phase transition. In this paper, we 
show that the gapped phases in isotropic two-leg spin lad- 
ders divide into two topologically distinct classes. This 
implies that phase transitions must necessarily occur if 
we try to interpolate between models belonging to differ- 
ent topological classes. Although not all gapped phases 
are directly equivalent to the Haldane phase of the spin-1 
chain, all possess some kind of string order. 

The rest of the paper is organized as follows. In Sec. II 
we introduce the spin ladder models that we will con- 
sider, and we recapitulate briefly what is known about 
these models. In Sec. HI we discuss the relationship be- 
tween valence bond states and string order, and also the 
possibility of phase transitions between spin models with 
topologically different string order. A bosonization treat- 



ment of the various models is presented in Sec. IV, and 
a discussion of the results obtained is given in Sec. V. 
Finally, Sec. VI gives a summary of our results and an 
outlook for future work. 



II. THE MODELS 

We begin with two antiferromagnetic spin-1/2 Heisen- 
berg chains with the Hamiltonian 

N 
'Ho = J 2_^ \^iA ' ■5*4+1,1 + Si^2 ■ '5'i+l,2 j , (2) 



where Si^i {Si^2) is the spin operator at site i on chain 1 
(chain 2). We will consider various forms for the inter- 
chain coupling. 

The interchain coupling 



N 



Til — Jl 2_^ ^iA ■ "^^2 



(3) 



describes the usual rung coupling along the legs of the 
ladder. This type of ladder is shown in Fig. g. (We will 
simply refer to this ladder model as a ladder.) 
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FIG. 2. The usual two-leg ladder. 

When Jl is strongly ferromagnetic ( Ji < and | Ji| ^ 
J) the two spins on the rung form a triplet, the singlet 
being much higher in energy. In this limit, the ladder 
behaves like a spin-1 chain, and hencCj-the spectrum is 
gapped. However, it has been shownO that a gap is 
generated by an arbitrarily small ferromagnetic coupling. 
Therefore, it appears that weak and strong ferromagnetic 
coupling are continuously related. 



FIG. 3. A typical configuration in the RVB state of the 
antiferromagnetic ladder. Solid lines represent singlet bonds. 

When Jl is strongly antiferromagnetic (Ji > and 
Jl » J), the ground state is essentially a product of 
rung singlets with a gap to magnon excitations. When 
Jl = J, it was shown that the ground state is well de- 
scribed by a nearest neighbor resonating valence bond 
(RVB) state and has a gap to the excited statesJlj A 
typical configuration of the RVB state is shown in Fig^ i 
Similar to the ferromagnetic case, it has been showrEjQ 
that the spectrum is gapped for arbitrarily small anti- 
ferromagnetic interchain coupling. Therefore, weak and 



strong antiferromagnetic coupling also seem to be con- 
tinuously related. 

We will also consider a ladder in which the interchain 
coupling is along plaquette diagonals 

N 

'H2 = J2 2_^ ^i,i ■ '^2+1:2 + Si^2 • (Si+ia . (4) 

j=l 

Together with Hq this gives the ladder shown in Fig. ^. 
(We will refer to this ladder as a diagonal ladder.) 



of the spins are bound into singlets and the ground state 
is unique. 
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FIG. 4. The diagonal ladder model. 

For J2 = J, this model is in fact the composite spin rep- 
resentation of a spin-1 chain. More precisely, by starting 
with the Hamiltonian of a spin-1 chain 



N 

1=1 



S*. 



'i+l 



and representing the spin-1 operator on site i as the sum 
of two spin-1/2 operators, Si = Si^i + Si^2, we find 

Ti = Ho + H2 , 

with J2 = J. In the composite spin representation, the 
total spin of each rung commutes with the Hamiltonian, 
so the eigenstates can all be classified by the total spin 
on each rung. The set of eigenstates with only triplets on 
all of the rungs corresponds to the spectrum of the spin- 
1 Heisenberg chain. Hence, the low energy spectrum of 
the composite |Sein representation is identical to that of 
a spin-1 chain.Ea The ground state of this model is well 
described by the AKLT statejj a typical configuration 
of whieii is shown in Fig. 0. Here again, it has been 
shownll2l that a gap appears for arbitrarily small J2 > 0- 
It has also been showrO that the spectrum is gapped for 
J2 < 0; in fact, the gap is generated for arbitrarily small 
J2 < 0. Therefore, it appears that weak and interme- 
diate coupling are continuously related for both J2 > 
and J2 < 0, however the gap vanishes at J2 = 0. 




FIG. 5. A typical configuration in the AKLT state of the 
diagonal ladder. Solid lines represent singlet bonds. 

As one can see from Fig. pi for a finite system there 
are effectively free spin-l/2's at the ends of the ladder, 
which are responsible for a four-fold degenerate ground 
state. When periodic boundary conditions are used, all 



FIG. 6. The zig-zag ladder shown as (a) a chain with first 
and second neighbor interactions (b) a ladder. 

Finally, we will consider an interchain coupling similar 
to 7^2 , but with only one of the diagonal couplings 



N 



H,3 — J3 2_^ 'S'i+l,! • Si^2- 



(5) 



i=l 



We consider such an interchain coupling because, for 
Ji = J3, we can write this model as a spin-1/2 Heisen- 
berg chain with first and second neighbor interactions 

"H = "Ho + Ti-i + 7Y3 

2N 2N 



(See Fig. 0). (We will refer to this type of ladder [for 
Ji — Js] as a zig-zag ladder.) 
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FIG. 7. The two degenerate ground states of the zig-zag 
ladder at the Majumdar- Ghosh point. Solid lines represent 
singlet bonds. 

The zig-zag ladder is knownt3 to have a quantum 
phase transition from critical, gapless behavior to a 
spontaneously dimerized gapped phase as J increases at 
J = ik^AUs. Moreover, at the special Majumdar-Ghosh 
pointtil (J = O.5J3), the ground state is known exactly. 
It is two-fold degenerate in the thermodynamic limit, and 
each of the ground states is a sequence of decoupled sin- 
glets, as shown in Fig. Q In fact, it has been shown that 
for 0.241 < J/ Ja^sj^xjo, the zig-zag ladder is gapped and 
has dimer order .E£ftj Therefore, it appears that the en- 
tire range 0.241 < J/J3 < 00 is continuously related to 
the Majumdar-Ghosh point. 



III. VALENCE BOND STATES AND STRING 
ORDER 

As previously mentioned, all of our models have short- 
ranged (SR) valence bond (VB) ground states, and all of 
our models exhibit string order. In this section, we will 
argue that there are two topologically distinct types of 
string order and that these two types of string order are 
intimately related to the VB structure. 

In general, any singlet state of an SU(2) symmetric 
model can be represented in terms of VB's. The SR-VB 
ground states of gapped spin liquids are, in the typical 
case, a linear combination of a large number of VB config- 
urations, in which the probability to find a longer ranged 
VB is exponentially small. The gap in the spectrum is 
related to the finite energy needed to break a VB. On the 
other hand, systems with a gapless spectrum necessarily 
contain longer VB's as well. 

In order to see the connection between the VB struc- 
ture and string order, we first consider the diagonal lad- 
der with ,Ji = J. The ground state is well described by 
the RVB picture of the AKLT state. A typical config- 
uration was shown in Fig. |5|. One particular spin con- 
figuration of Fig. H is shown in Fig. H(a). Suppose we 
add the z component of the spins on the same rung, as 
shown in Fig. @(b). The total Sf can take on the val- 
ues 1,0,-1. Considering this sequence, if we remove all 
sites with Sf = 0, the remaining sites have Neel order 
(i.e. there is string order). If, on the other hand, the 
z components of the spins along plaquette diagonals are 
added, as shown in Fig. ||(c), there is no string order. 
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FIG. 8. (a) One particular spin configuration of Fig. ^ 
(b) Dashed lines encircling rungs: S't — Si^i + S,?2- Notice 
that after we remove all sites with Sf = 0, the remaining 
sites have Neel order, (c) Dashed lines encircling diagonals: 
S! = Sf+i,! + Sl2- (Figure adopted from Ref. |.) 

Now consider the antiferromagnetic ladder with Ji — 
J. The ground state is well described by a nearest- 
neighbor RVB state, for which a typical configuration 
was shown in Fig. 0. One particular spin configuration 
is shown in Fig. |9|(a). Suppose we look at the z compo- 
nent of the total spin on a rung, as shown in Fig. P(b); 
the state has no string order. If, however, we consider 



the z component of the total spin along plaquette diag- 
onals, as shown in Kg. |^(c), string order is found. It 
was shown by WhiteB that there is a 96.2% probability 
of finding triplets akjag plaquette diagonals, and it was 
verified numericalljoEl that these, triplets in fact exhibit 
string order. It was also showrocl that no string order is 
found if the total spin along rungs is considered. 
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FIG. 9. (a) One particular spin configuration of Fig. g. (b) 
Dashed lines encircling rungs: Sf — Sf^i + Sf^- (c) Dashed 
lines encircling plaquette diagonals: Sf = 'S'f+i,! -I- Si2- No- 
tice that after we remove all sites with S^ = 0, the remaining 
sites have Neel order. (Figure adopted from Ref. |9|.) 



As mentioned before, the ladder with ferromagnetic 
coupling along the rungs is continuously related to the 
true spin-1 chain. Therefore, it has an AKLT-like ground 
state and string order due to triplets along the rungs. On 
the other hand, the ground state of the diagonal ladder 
with ferromagnetic interchain coupling has an RVB-like 
ground state, similar to that of the ladder with antiferro- 
magnetic interchain coupling; as discussed above string 
order is due to triplets along plaquette diagonals. 

Finally, consider the zig-zag ladder at the Majumdar- 
Ghosh point. The exact ground state is given by a 
product of decoupled singlets, as shown in Fig. 0. One 
particular spin configuration for Fig. OTa) is shown in 
Fig. |l^(a) . Adding the z components of the spins on the 
same rung, as shown in Fig. nn(b), it is obvious that we 
always get zero. However, if we add the spins along pla- 
quette diagonals, as shown in Fig. |l^(c), we find string 
order. Now, however, consider the ground state shown 
in Fig. M(b). One particular spin configuration is shown 
in Fig. |ll|(a). Then, as shown in Fig. |l|(b), if we add 
the spins along the rungs, we find string order. However, 
if we add the spins along plaquette diagonals, as shown 
in Fig. |ll|(c), we always get zero. Jjhis result is not sur- 
prising, since it was already shownHlthat the Majumdar- 
Ghosh ground state has perfect string order. However, 
it appears that, depending on which of the two degen- 
erate ground states is actually realized, string order can 
be due to the spins along the rungs or to the spins along 
plaquette diagonals. 
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FIG. 10. (a) One particular spin configuration of Fig. |7|(a). 
(b) Dashed lines encircling rungs: SI = SI i + Sf 2- (c) Dashed 
lines encircling plaquette diagonals: Sf = 'S'i+i,! + Si^2- No- 
tice that after we remove all sites with SI = 0, the remaining 
sites have Neel order. 



FIG. 11. (a) One particular spin configuration of Fig. M(b). 
(b) Dashed lines encircling rungs; Si = Si^ + S'f_2- Notice 
that after we remove all sites with SI — 0, the remaining 
sites have Neel order, (c) Dashed lines encircling plaquette 
diagonals: Sf — Sf+i.i + 5'f_2- 



Motivated by the above examples, similar to Rcfs. and 0, we introduce two string order parameters 

a"dd = - i^j™^ /(^M + 3^2) exp u J2 isr,i + sr,2y\ isf.i + sf,^)\ , 

(^"ven = - i^ijm^ /(5f+i,i + Sr,2) exp ^TT ^ (8^+,,, + 5^) j (^jVi.i + Sf,2)j ■ 



(6) 



(The names Oodd and Ocvcn will be made clear below.) We saw that when one of the order parameters is finite, the 
other vanishes. So, the antiferromagnetic ladder has Oovcn 7^ 0; the composite spin model has Oodd 7^ 0; the zig-zag 
ladder at the Majumdar-Ghosh point is special - it can have either Oovon 7^ or Oodd 7^ 0, depending on which of 
the two degenerate ground states is actually realized. However, it cannot have both finite simultaneously. 

Let us now consider the topology of the VB's in the above examples; an interesting pattern emerges. 
If we count the number of VB's crossing an arbitrary vertical line, we find that this number is always 
even for the ground state configurations of the antiferromagnetic ladder (Fig. 0), while it is always odd 
for the diagonal ladder at J2 = J (Fig. ||). For the zig-zag ladder at the Majumdar-Ghosh point, 
the number of VB's crossing an arbitrary vertical line depends on which of the two degenerate ground 
states we consider: it is even for the state in Fig. 0(a), while it is odd for the state in Fig. 0(b). 
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FIG. 12. Solid lines connecting sites represent VB's. (a) 
Qy — odd: the number of VB's crossing a vertical line is odd. 
(b) Qy — even: the number of VB's crossing a vertical line is 
even. 



These examples are special cases of a more genepJ clas- 
sification of SR-VB states. It has been showrJSI that 



for SR-VB configurations on a two-dimensional square 
lattice, two topological numbers, Qx = even/ odd and 
Qy = even/odd, can be defined. They are determined by 
the parity of the number of SR-VB 's crossing arbitrary 
horizontal and vertical lines parallel to the x- and y-axes, 
respectively. In the case of 2-leg ladders, only Qy is rel- 
evant; Qy can be either even or odd, as illustrated more 
generally in Fig. |l^. Hence, there is a topological number 
which distinguishes between whether the number of SR- 
VB's cut by a vertical line is even {Qy — even) or odd 
{Qy = odd). For any finite size system, the even and odd 
sectors are coupled as long as there are VB's with length 
comparable to the system size. However, when the sys- 
tem is gapped and thus has a SR-VB ground state, the 
tunneling amplitude between the two sectors goes to zero 
exponentially fast as L ^ 00, and the ground state is a 
pure Qy = even or Qy = odd state in the thermodynamic 
limit. Note that in long-ranged VB ground states of gap- 
less models, the even and odd sectors remain coupled; 



hence, no such topological distinction is possible. 

It is also worth noting that for open boundary condi- 
tions Qy = odd ground states have spin-l/2's locaHzed at 
the ends of the ladder, while Qy = even states do not. As 
is obvious from Fig. |l2|, these end spins occur for topo- 
logical reasons. Their presence or absence is probably 
the simplest way to determine Qy. 

From the above examples, the (topological) parity of 
the SR-VB ground state and the type of string order seem 
to be intimately related. It appears that ground states 
with Qy = odd have Oodd string order, while ground 
states with Qy = even have Oeven string order. 

Now suppose we smoothly vary the parameters of the 
Hamiltonian, such that we interpolate between models 
belonging to different topological classes; a phase transi- 
tion necessarily occurs. A priori this transition could be 
either first order or second order, depending on the actual 
path in parameter space. When the transition is second 
order, the string order parameters vanish at the transi- 
tion point and the ground state becomes a long-ranged 
VB state. In the next section we analyze this problem in 
the weak coupling limit using bosonization. 



IV. WEAK COUPLING ANALYSIS: 
BOSONIZATION 

In the bosonization treatment of the ladder model 
shown in Fig. ^, we start with two decoupled spin-1/2 
chains and treat the interchain couplings perturbatively. 
Our conventions, as well as the bosonization of spin lad- 
ders, are presented in detail in Ref. p^. 




FIG. 13. Most general ladder model which we will consider, 
Ti. = Ho + T~ii + T~i2 + Ha- 

The isotropic spin-1/2 Heisenberg chain is known to be 
critical. The effective Hamiltonian for long wavelength 
excitations is 



Us 



= 1/2 



= — I dx 



KTf + ^{d^<^f 



(7) 



where the bosonic phase field, $, and its conjugate mo- 
mentum, n, satisfy the commutation relation 



mx),n{y)]^^5{x-y)- 



(8) 



For an isotropic antiferromagnetic spin chain, K — 1/2. 
We will also need the bosonized form of the spin opera- 
tors. They areE3 



^+ ^ S^ ^ exp (-^^/^9) 



'i{-KX / a) 



-I- cos(v47r$) 



a 



S^{x) =^-^= 1 9,$ + e^(Wa)!WV4-£) 



(9) 



where the dual field, O, is related to Ft by 11 = 9x6. 

For the ladder, we simply attach a chain index to our 
fields. Therefore, 



Ti-a ^ ^ I dx 
+ ^ I dx 



ifn2 + l(9,$i)' 
Kul + -i(a,<i>2)^ 



(10) 



To bosonize the interchain coupling, we write the spin 
operators in terms of uniform and staggered components 
as 



S,(x-)=J»(a;) + (-l)"/"n,(a;). 



(11) 



We find 

Hi ^ Ji I dx [Ji(x) • 32{x) +ni(a;) • n2(a;)] , 
7^2 ^ 2J2 / dx [3i{x) ■ 32{x) ~ Ta.i{x) ■ n2{x)] , 
Hs-Jaf dx [Ji(x) • J2(x) - ni(a;) • n2(:E)] . (12) 

Inserting the expressions in Eq. (M into Eq. (O) gives 

— ^ \gl cos (V4^($i H 
(27ra)^ L V 

-H 52 cos f V47r($i -$2) 

4 
J] 



n^ 



5^cos(v^(ei-e2))] 



dxdx^idx^2 



dx 



{2nay 



[glcos{V^{ei-Q2)) 



X cos (\/47r(<i>i + $2 

+ 5^cos(0F(ei-e2)) 

X cos (\/47r(<i>i - $2 



(13) 



The gl^ gl, and g^ terms come from ni(x) • n2(a;), while 
the J^, g\, and gl terms come from Ji(a;) • 32{x). 
For Hi: 



gl = -2Ji; gl = 2Ji; gl = 27rJi, 
J^ = Ji; gl = ttJi; gl = ttJi . 



For H2: 



gl^AJ2\ gi^-iJi; gj ^ -4ttJi, 
Ji=2J2; gl^2nJ2; ffi = Stt J2 . 



(14) 



(15) 



For H3: 



gl = 2J3; gl = -2J3; g| = -27rJ3, 

^i = -/a; ffl = ttJs; .95 = TrJg . (16) 

It is useful to define the fields 

$,,, = -i= ($1 ± $2) , e,,a = -^ (61 ± 62) . (17) 

In terms of these fields our Hamiltonian, H = Ho +7ii + 
7^2+^3, is 



[cos (V^^Qa) cos {VS^<^a)] = '2Ka + l/{2Ka). There- 
fore, gi will grow at large distances for Kg < 1; 32 
will grow for Ka < 1; <?3 will grow for Ka > 1/4; 
(74 will grow for 2Ks + l/(2Ka) < 2; 55 will grow for 

2Ka + l/{2Ka) < 2. 

In what follows, we will consider the phases and tran- 
sitions that occur when we vary Ji, J2, and J3. In or- 
der to make things more tractable, we will consider two- 
dimensional slices in the full Ji — J2 — J3 space. 
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dx 



(2TTay 



■ cos 



(V8^$s) 



dx 






+54 



dx 
(27ra)2 
dx 



52 cos 



+55 



(27ra)2 

dx 
(2^ 



(V8^$a) + 53 cos (V2^ea) 
('\/2^ea) cos (V8^<i>s] 

(^V2^ea) cos (V8^$a) , (18) 



where 



51 = -2 Ji + 4J2 + 2 J3 , 32 = 2 Ji - 4J2 - 2 J3 , 

53 = 27r Ji - 47r J2 - 27r J3 , 54 = ttJi + 27r J2 -t- tt J3 , 

55 = TT Ji + 27r J2 + TT J3 , 

Jj_ = Ji + 2 J2 + J3 . (19) 



Also, 



K.=K I 



K„^K I 



KJ, 



KJ_i_ 



-1/2 



-1/2 



u 1-K 



i^J, 



u 1 



KJ_i_ 



1/2 

(20) 
1/2 



For Jj^ <C 1 we have 

KJ±\ 



K^^K 1- 



2uTr 



■ 



2UTX J 



2UTT I 



U 1 



2ut: 



(21) 



We are interested in whether or not the inter- 
chain coupling causes a gap in the excitation spec- 
trum. Therefore, we would like to identify the rel- 
evant operators; these operators will "pin" their ar- 
guments, thus causing gaps to appear. To do this 
we consider the scaling dimensions of the operators 
in the interchain couplingB3'E3 The scaling dimensions 
of the operators are the following: [cos (V87r$s)] = 

2Ks; [cos(V8^$a)] = 2Ka; [cos{V2^ea)] = 
l/{2Ka); [cos {V2^&a) cos (\/8^$.)] = 2Ks + l/{2Ka); 
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FIG. 14. Phase diagram for Ji / and J2 / with J3 = 0. 



A. Ji / and J2 7^ with J3 = 

In this case, for J2 = and Ji = J we recover the usual 
antiferromagnetic ladder; for J2 = and Ji ^- — 00 we 
recover a spin-1 chain, where the spins on each rung form 
an S* = 1. Similarly, for Ji = and J2 = J we recover 
the composite spin representation for a spin-1 chain, and 
it was previously shown that the composite spin repre- 
sentation haa-|the same low energy physics as the true 
spin-1 chain.113 

The phase diagram in the J\ — Ji plane is shown 
in Fig. n3. In region / the g\ and (?3 terms are the 
most relevant. Therefore, <I>s and 0a are pinned with 
($s) = (2n + l)v^7;78 and (9^) = toV2^. Without loss 
of generality, we can choose n = and m = 0. This 
gives ($s) = \/7r/8 and (6a) = 0. In region II the g\ 
and g-i terms are again the most relevant. Again, $«, 
and Oa are pinned. Now, however, ($s) = n^JiT /2 and 
(0a) = (2m + l)y^r/2. Choosing n = and m — 0, 
($s) = and (6a) = ^7r/2. In region /// the gi and gi 
terms are the most rele vant . Therefore, $s and <I> a are 
pinned with ($s) = n^7r/2 and (^a) = (2m -I- 1)^/77/8. 
Choosing n = and m = gives {^s) — and 
($a) = -\/7r/8. In region /]/, similar to region ///, the 
gi and gi terms are the most relevant, so $<, and ^a are 
pinned. However, in this region (<I>s) — {2n + l)-\/7r/8 
and ($a) = m-\/7r/2. Choosing n = and m. = 0, 
($s) = ^tt/S and ($a) = 0. There are also two special 



lines in the phase diagram. Along the line J2 = J'l/S, the 
111 •n2 terms vanish and only the Ji • J2 terms remain. For 
Ji, J2 < 0, the system is gapless; for Ji, J2 > 0, the 34 
term is marginally relevant and the system is gapped. 
However, the ground state is two-fold degenerate (for 
Ji,J2 > 0): ($,) = (271 + 1)V^, (Qa) = 2mV^ 
or ($5) = 1n.jTx]%, (Ga) = (2m + 1)^7^72- Choosing 
n = and m = 0, we have ($s) = a/tt/S, (9a) = 
or ($s) = 0, (0a) = Y^7r/2. The other special line is 
J2 = — Ji/2. Along this line the Ji • J2 terms vanish and 
the ni • n2 terms all have the same scaling dimension = 
1. Along this line the spectrum is gapped. This line will 
be discussed in greater detail in Sec. V. 



h=\ 




T he p hase diagram in the J\ — J3 plane is shown in 
Fig. 15. The regions /, //, ///, and IV have properties 
identical to the J\ — J2 phase diagram discussed above. 
The line J3 = —J\ has properties identical to the line 
Ji — —Ji/2 discussed above, and will be discussed in 
greater detail in Sec. V. The line, Ji = J3 is special. As 
pointed out by Nersesyan et al.fA we must be careful of 
the "twist" operators which appear. Along this line, the 
interchain coupling of the staggered components can be 
written as 



Explicitly, the terms are 

/(Jlx / 

— — -yaj;$asin V87r$ 

+ 52 



(22) 



53 



dx 

(27ra)2 

dx 
(27ra)2 



9j,$sSinf V87r$a 
aj;e,5sinfV27reaj , 



(23) 



FIG. 15. Phase diagram for Ji / and J3 / with J2 = 0. 



where g'^, g'l, g'^ '^ J3. These terms are subtle because 
they have non-zero conformal spin. As pointed out in 
Ref. |2^, a seemingly irrelevant operator with non-zero 
conformal spin can generate relevant operators. However, 
since we are only considering the SU(2) symmetric case, 
the terms generated are less relevant than the Ji-J2 terms 
already present. Therefore, similar to the line J2 = Ji/2, 
for J3 < the system is gapless; for J3 > the g^ term 
is marginally relevant and the spectrum is gapped, with 
the ground state being two-fold degenerate. 



V. DISCUSSION OF THE RESULTS 



B. Ji / and J3 / with J2 = 

In this case, for J3 = and Ji = J we recover the usual 
antiferromagnetic ladder; for J3 = and Ji — > —00 we 
recover a true spin-1 chain. For Ji = J3 7^ we have a 
zig-zag ladder; in particular, for Ji = J3 = 2J we have 
the Majumdar- Ghosh point where the ground state is 
dimerized with two-fold degeneracy. 



As discussed in Sec. HI, the two-leg ladder models we 
have considered can have either O'^^^ or Ocvon string 
order, but not both simultaneously. In our case, since 
we are only considering SU(2) symmetric models, O^ = 
Qy — Qz^ Therefore, for simplicity, we focus only on O^. 

To bosonize the string order parameter, we first write 
it in a more convenient form. Using the identity 
exp (iirS") = 2iS^ , we can write 



o- 



dd 



lini 



4 \t-j\- 



1 



Ceven = 7 lip 

4 l«-j|^oo 



'e^'^S-i _^g»7rS, 



Bosonizing O^^^ and Oovcn gives 



Of 



lim ( e 



e"'^'+i'i 



exp Ln J2 (^U + ^"2) ] (e*"^-^ + e^"^-^) ) , 
e*-^°^) exp Ltt J2 (^r+1,1 + ^"2) ) (e"^?+M +e-^°^) 



(24) 



l=i+l 



iV27r*s (k) -i\/27r*s (j/) \ 



(25) 



• cvcn/odd , '^^Y^ \ "^ ■-- / 

We see that all we need is for $s to get pinned to have 



string order. The operators for both O^^^ and O^^^^-^ have 
the same bosonized form because the nonlocal string op- 
erator makes the continuum limit insensitive to physics 
occurring on the order of a single lattice spacing, such as 
whether triplets lie predominantly along rungs or along 



diagonals. Therefore, the bosonized string order param- 
eter tells us that we have string order, but it does not 
tell us in which topological sector the order exists. Tak- 
ing into account the physical picture we get from the 



VB states in Sec. Ill, we can understand the various re- 
gions and transition lines which have been obtained in the 
phase diagrams by bosonization. The results are summa- 
rized in Tabic ||. 



TABLE I. Results for the Ji — J2 and Ji — J3 phase diagrams. 
/ II III 



IV 



<^s 




($,) = ^it/8 ($,) = 




($,) - ($,) - VV8 


<^a,Qa 




i&a) = (6,) = v^ 




($a> = VV8 ($a> - 


Order Parameter 




^odd ^cvcn 




^evcn ^odd 


J2 = Ji/2: 




J2,J3 > — first order transition 




J2,J3<0 — second order transition 


(J3 = Ji) 


(*.) = 


Vtt/S , {Qa) = or ($,) = , (e„> = , 


A/2 


$3 and $a critical 


J2 = -Ji/2: 




level crossing in the excited states 






[Js = ~Ji) 











A. Ji / and J2 / with J3 = 

The line Ji = with J2 > is continuously related to 
the composite spin model; the composite spin model has 
Oodd string order. Therefore, it appears that region / is 
continuously related to this model, and hence has Oodd 
string order. The line J2 = with Ji > is continu- 
ously related to the usual antiferromagnetic ladder; the 
antiferromagnetic ladder has Oovcn string order. There- 
fore, it appears that region // is continuously related to 
this model, and hence has Ocvcn string order. Along the 
line Ji ~ with J2 < 0, we have ferromagnetic inter- 
chain coupling along plaquette diagonals. For | J2I ~ Ji 
the ground state is similar to the RVB state of the lad- 
der with antiferromagnetic interchain coupling. Hence, 
this model has Ocvcn string order. It appears that region 
/// is continuously related to this model, and hence has 
Oeven String Order. Finally, the line J2 — with Ji < 
is continuously related to the spin-1 chain in which the 
spins on a rung form an effective S = \. Since the ground 
state of the spin-1 chain is described by the AKLT state, 
this model has Oodd string order. Therefore, it appears 
that region IV is continuously related to this model, and 
hence has Oodd string order. We see that a transition 
between Oodd and Ocvcn string order occurs along the 
line J2 = Ji/2. For J2 < the transition is second order; 
for J2 > 0, there is a marginally relevant operator which 
drives the transition first order. The line J2 = — Ji/2 is 
interesting, so we discuss it in detail. 

Along the line J2 — ~J\I2 there is a change in the 
properties of the system: above the line, 0a is pinned; 
below the line, $a is pinned. However, we believe this 
is a level crossing in the excited states; the properties of 
the ground state remain the same. Hence, the system 
does not undergo a phase transition when we cross this 



line. To show this, it is useful to express Ha in terms of 
Majorana fermions. We begin on the line J 2 = — Ji/2; 
along this line, the Ji • J2 terms vanish and the rii • r\.i 
terms all have the same scaling dimension. Rescaling our 
fields. 



(26) 



Ha has the form 

H„ = ^|dx[H2-i-(5A)2] 

[y^^^ + g-i cos (Vl^Oa 



52 COS 



(27ra)2 
Using thato 

= -iu \ dx ^jj5^-0fl, - tpl^d,j:ip_ 
— cos(\/47r$) = i (tp]itpL - -^iV'fl) 7 

j- cos(V4^e) = I (vljVl - i^Wn) , 

Ha can be written as 

Ha = -iUa / dx Njljid^1pa,R " i^l^L^^x^a 
92 [lplR'^a,L - i'i.Li'a 



(27) 



(28) 



dx 
2'na 



'93 



(29) 



Now introduce two independent Majorana fermions, ^ 
and 77, defined by 

i^a,R = —f={^B. + IVb) , tpa,L ^ —^{^L + IVl) ■ (30) 

Finally, Ha can be written as 



+ iVRdxTiR - IlLdxVL)] 
+ 2 / ^ [(52 - 93)^R^L + (52 + 93)vrVl] 



(31) 



This is the Hamiltonian for two massive Majorana 
fermions. As is well known, massive Majorana fermions 
describe the long distance properties of the Ising model 
away from criticality. 

For J2 w -Ji/2 (i.e. J2 = -Ji/2 + S, \6\ < 1), the 
Ji • J2 terms do not vanish. However, very close to the 
line J2 = — >/i/2, we can still write Ha in terms of Ma- 
jorana fermions. The Ji • J2 terms can be written as 
four-fermion interactions which just renormalize the ve- 
locity and fermion massesE3. The key thing to notice is 
that when we cross the line J2 = — Ji/2, the values of the 
fermion masses change, but their signs do not change. It 
is well known that the (Majorana) fermion mass chang- 
ing sign corresponds to the order-disorder transition of 
the Ising model. Since there is no change in sign when 
we cross the line J2 = — Ji/2, the structure of the ground 
state does not appear to change. Therefore, we interpret 
the change from Qa being pinned to ^a being pinned as 
a level crossing in the excited states. Hence, the sys- 
tem does not appear to undergo a phase transition when 
crossing this line. 



B. Ji / and J3 / with J2 = 

The line Ji > with J3 = is continuously related 
to the usual antiferromagnetic ladder; the antiferromag- 
netic ladder has Oevcn string order. Therefore, it appears 
that region // is continuously related to this model, and 
hence has Oovcn string order. The line Ji < with J3 — 
is continuously related to the spin-1 chain in which the 
spins on a rung form an effective S = I. Since the ground 
state of the spin-1 chain is described by the AKLT state, 
this model has Oodd string order. Therefore it appears 
that region IV is continuously related to this model, and 
hence has Oodd string order. Along the line Ji — 0, the 
coupling along the rungs is zero and only the diagonal 
interchain coupling, J3, is nonzero. This is similar to the 
case when only Ji ^ 0, except with chain-1 shifted to 
the right by one lattice spacing. Therefore, J3 > is 
similar to the usual antiferromagnetic ladder and J3 < 
is similar to the spin-1 chain, except with chain-1 shifted 
to the right by one lattice spacing. Hence, J3 > has 
Codd string order, and J3 < has Ocvcn string order. It 
appears that region / is continuously related to the line 
J3 > and that region /// is continuously related to the 



line J3 < 0. Therefore, region / has Codd string order 
and region III has Ocvon string order. A phase transition 
occurs along the line J3 — Ji. For J3 < the transition is 
second order, while for J3 > a marginally relevant op- 
erator drives the transition first order. Similar to the line 
J2 = — Ji/2, the system changes character when crossing 
the line ,73 = — Ji. Above the line, Oq is pinned; below 
the line <f>a, is pinned. Similarly, we believe that there 
is no phase transition as we cross this line; it is a level 
crossing in the excited states. 

It is interesting to note that in our model, the zig-zag 
ladder (i.e. the line J3 = Ji) is actually a transition 
line. For J3 > the line is a first order transition line, 
while for J3 < the transition is second order. There- 
fore, the Majumdar- Ghosh point actually lies on a first 
order transition line. 



VI. CONCLUDING REMARKS 

In this paper we studied the gapped phases in two-leg 
spin ladders. The ground states of these ladders are well 
described by SR-VB states. There are two topologically 
distinct classes characterized by whether the number of 
VB's cut by a vertical line is even {Qy — even) or odd 
{Qy — odd). Note that this classification of Qy — even 
and Qy = odd can be used for even-leg ladders but not 
for odd-leg ladders. For odd-leg ladders, one gets an 
even-odd alternation, as shown schematically in Fig. Hq. 
This even-odd alternation implies a two-fold degener- 
ate ground state, consistent with the Lieb-Schultz-Mattis 
theorem£d 




odd 



FIG. 16. A VB configuration which could occur for a 
three-leg ladder, provided the couplings are chosen so that 
the ground state is described by a SR-VB state. Notice that 
the number of bonds crossing a vertical line alternates. 

Associated with Qy = even and Qy = odd, we con- 
sidered the "even" and "odd" string order parameters 
of Eq. g for the ladder model shown in Fig. |l^. Using 
known results for particular values of the coupling con- 
stants along with bosonization, we obtained the phase 
diagrams in the Ji — J2 and Ji — J3 planes, shown in 
Figs. 04 and la, respectively. While these results cover 
only parts of the Ji — J2 — J3 parameter space, we believe 
that the association of the string order parameters C"vcn 
and ©odd with Qy = even and Qy = odd is appropriate 
for this model in general. 

We should emphasize that this classification of Qy = 
even and Qy = odd relies on the possibility of writing the 
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'^^' and spin ladders with 4-spin plaquette 
both having a non-Haldane-like dimerized 



singlet ground state of the ladder as a superposition of 
VB configurations. This may break down in anisotropic 
models, and so by introducing anisotropic couplings the 
two topological sectors may be coupled. 

Another interesting problem is related to the 
spin-1 chaip^ with bilinear and biquadratic exchange 
interactio 
couplings Jl; 

phase. Although Qy is still a good topological number 
for the dimerized ground state, it is not clear if string or- 
der is simply due to the short-ranged nature of the VB's 
and survives the transition from the Haldane phase to 
the dimerized phase. 

It is clear from this analysis that the apparently fea- 
tureless spin liquid phase of spin-gapped two-leg ladders 
actually has a rich underlying topological structure. It 
remains to be seen what role these ideas may play in 
the doped systems. More precisely, does this topologi- 
cal structure survive when the system is doped, and is 
pairing ultimately related to the topological structure? 
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